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Проанализированьт основньюе математическиє модели ограниченной задами логарифмического потенциала 
для контактной поверхности. Отмечено, что параметрическая зависимость решений задачи от плотности актив- 
ного слоя и от "средней глубиньї" или от асимптотьт контакта являєтся следствиєм изучаємьтх моделей сложной 
средь, а не принципиальной характеризацией задачи для неограниченньх тел в отличие от задачи для ограничен- 
ньїх тел. Детально исследовано нелинейное интегральное уравнение для контактной поверхности в комплексной 
плоскости, которое зависит только от одного параметра - плотности гравитирующего слоя. Установленьт условия 
корректной разрешимости уравнения. Доказаньт локальньте теоремьт единственности, существования и устойчи- 
вости его решения на определенном компактном множестве. 
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Вступ. Обернена задача потенціалу для контактної поверхні сформульовані давно 1 її досліджен- 
ня висвітлені в численній геофізичній літературі, яка час від часу поновлюється. Разом з тим і досі не 
з'ясовано питання, з якими особливостями задачі пов'язано розмаїття її описів: чи то з неозорою мно- 
жиною моделей задачі; чи то з наявністю теоретичних проблем її розв'язності, що важко піддаються до- 
ність у літературі детального аналізу постановок задачі 1 методів її розв'язування, з якого випливала б 
відповідь на поставлене питання, спонукає нас до спроби з'ясувати ситуацію, що склалася. 

Проте ми не ставимо перед собою мети надати вичерпний огляд всієї літератури, в якій з тою чи 
іншою мірою повноти й строгості обмірковано проблеми коректної розв'язності й чисельної реалізації 
розв'язків задачі для контактної поверхні і вивчено ефективність їх використання для тлумачення при- 
роди гравітаційних та магнітних аномалій. Ми обмежились якісним аналізом основних, на нашу думку, 
праць, що проливають світло на стан і рівень вивченості математичної моделі задачі, залишаючи поза 
увагою багато цікавих питань прикладного, історичного, можливо, й пріоритетного плану. 

Основні моделі задачі. Задачу про відновлення межі поділу двох однорідних середовищ уперше 
сформулював 1 наближено розв'язав Б.В. Нумеровим в 1930 р. у такій постановці |1). Відомо, що під 
плоскою поверхнею Землі ОД: у «0 є одна кон- 
тактна границя ОР, : б з СО»), -о«хаю,що 
" 


розділяє однорідні шари легких "верхніх " порід з 
густиною С; - соп5і від важчих "нижніх" порід з 
густиною с, «соп5і. Окрім різниці с -«с0,-6» 
густин порід задано значення сили тяжіння 
0, (0), -о«х«со, і віддаль Й 202) між по- 


верхнями ОР, і ОР, у певній точці х , яку відра- 





ховують за нормаллю до ОР), (рис. 1). За цими ви- 


хідними даними потрібно визначити "поверхню" ОР, припустивши, що невідома функція б «С (х), 


яка описує межу поділу, є однозначною й обмеженою. Для визначення контакту в цих умовах Б. В. Ну- 
меров вивів нелінійне інтегральне рівняння 


їз Чеоорее |о) 0 (0) 
що | п о абЕ - ве 1 
сно пт оц зу й 
що залежить від двох параметрів - стрибка с густини на контакті 0 (х) і глибини й залягання контакту 


а 5 зл г ж є й 
в деякій довільно обраній точці х , яка не збігається, взагалі кажучи, з точкою х спостереження поля; / 
- гравітаційна стала. 


Наближений розв'язок цього рівняння Б. В. Нумеров отримав у вигляді 
сн (») піди Ане(х) 
за додаткової умови, що рельєф контактної границі не надто контрастний, в крайньому випадку такий, 
що коливання функції ДЛ, 0 (х) менше від /. За цієї умови однозначну функцію Дні (х) визначено як об- 
межений розв'язок лінійного інтегрального рівняння першого роду 


15 й М й 
иа 2 


яке залежить від двох параметрів с 1 й, тобто Ане(х ) З АС гой ) 
М.Р. Малкін |2, 3) задачу визначення контактної поверхні 0/,, яка розділяє два однорідні косо- 





намагнічені шари, за заданим магнітним потенціалом і стрибком ./ - (г ,) вектора намагніченості по- 


рід зводить до розв'язування нелінійного інтегрального рівняння, яке за умови вертикальної намагніче- 
ності порід цілком збігається з рівнянням (1). Наближений його розв'язок су, (х) п АНУ (х) одержано 


методом Б.В. Нумерова. Тому він залежить від трьох параметрів: двох природних - складових ./, і ./, 


стрибка вектора намагніченості порід 1 одного "привнесеного" - товщини Й виділеного однорідного ма- 
гнітного шару, тобто Дб (х) чай (о У ЯР. | ;)). 


у? 
Г.Рейнбой (4), порівнюючи лінеаризоване рівняння (2) для контактної поверхні з інтегралом 
Пуассона для подання гармонійної функції у(Х, у) в півплощині за значенням її на рівні у - Й 


17 ьо, 
й І. сохузі (гл, 


виявив, що лінеаризована задача для контактної поверхні еквівалентна задачі визначення гармонічної 
функції у(х, у) за значеннями у(20), тобто що Да (х) о УС п). Тим самим він фактично показав, що 





точний розв'язок інтегрального рівняння (2) повинен бути слідом гармонічної функції (цей факт ще 
довго залишався неусвідомленим; 1, окрім того, відкрив перспективний напрямок в теорії інтерпретації 
аномалій, пов'язаний з аналітичним продовженням полів у півпростір зосередження їх джерел (не тіль- 
ки необмежених типу шару, але й обмежених) (5-10). 

О. О. Заморєв ЦП, 12| дослідив в комплекссній множині 2 - х -іу задачу визначення контактної 
поверхні ОД), у вигляді скінченної ундуляції у(х) шля п), у(а) - У(Б) -Л, на скінченному інтервалі 
(а, Б| горизонтальної" прямої УС) й асоп5ії» 0, хє а, 5) (рис. 2). Проблему відновлення форми збу- 
рюючого тіла (виповненого масами однорідної гус- 


тини с 1 обмеженого, на відміну від нумерівського 
шару, прямою уєП і кривою 


уСх) -йз п), хє а, 5, па) -- ть) -0), тобто 
задачу визначення параметрів а,Б,й однозначної на 
інтервалі |а,»| обмеженої функції лох) за заданими 
на осі у - 0 значеннями (відповідним чином прока- 





ліброваними) похідної логарифмічного потенціалу 
С 6 5). О.О. Заморев звів до нелінійного рівняння, 


Рис. 2. Модель Заморева. 


2 яке можна записати в двох еквівалентних один од- 


Ь 5 
б-хчій 1-й (2 ) 1 
І шо ; 3 
вени Оти НЕ йно) баз ща) У 
зда С) суч 93) аг! 90, у) :9И(2,у) 
ліс да 2 дх бу 
Обидва ці подання наведено у статті |111, до того ж друге, яке детально вивчено, одержано з пер- 
шого інтегруванням частинами за умови, що та) са ть) - 0. В результаті аналізу другого подання нап- 


ному зображеннях: 











руженості ц(2) поля О. О. Заморєв знайшов його точне обернення. На жаль, ефективність цього обер- 
нення при розв'язуванні прикладних задач не перевірено ще й досі. Звернемо тепер увагу на перше зоб- 
раження напруженості поля. Неважко виявити, що дійсна її складова 


А Р я а х)у-Ейчнх), х є|а 
те се о о а) 


3 
за зовнішнім виглядом нагадує рівняння Б. В. Нумерова (1) без доданку / в правій частині, що відтво- 
рює за вибраної калібровки поля притягання однорідного шару товщиною /?20. Більше того, величина 
інтегралу в лівій частині рівняння (4) не змінюється на заморєвських тілах у разі заміни скінченних гра- 
ниць а 1 на нескінченні -оо і со, тому ліві частини рівнянь (1) і (4) за зовнішнім виглядом повністю збі- 
гаються, хоча і зображають в математичній символіці різні за суттю предметні моделі: перше - притя- 
гання необмеженим шаром між границями у- ПЛ і у-а (х), -о0« х «00; друге - притягання скінчен- 
ним тілом, обмеженим кривими у-п/ і укич лох) у Хе а, 5) ; та) - ть) - 0. Можливо, зовнішня по- 
дібність лівих частин рівнянь Б.В. Нумерова й О. О. Заморєва відігравала свою роль у подальших прик- 


ладних дослідженнях, серед яких осібно вирізняється праця |131. 
У цій праці, написаній в 1941 р. і виданій в 1947 р., О.О. Шванк помітив, що коли в рівнянні Б.В. 


Нумерова замість параметра Й підставити значення контакту С (х) в точці х спостереження поля, то 
з рівняння (1) виходить нове нелінійне інтегральне рівняння 


а В А, з 
З ей гир оон (5) 


яке за не дуже великих коливань шуканої функції 0 (х) можна розв'язувати методом послідовних наб- 





лижень. На одному окремому прикладі він показав, що обчислені "вручну" послідовні наближення 
я й (х)) контакту 0 (х) з першого 2, (х) - и(х) і до 6-го включно мають тенденцію до збіжності. Проте 


зауваження О.О. Шванка, засноване на виводі Б.В. Нумерова, залишилось, на жаль, без уваги, можливо, 
через те, що і в самого автора пропозиції не склалося чіткої думки про єдино можливий функціональний 
зв'язок поміж контактом і полем. 

Дійсно, в 5 2 розд. 85 (стор. 361-374) його книги |13| для поля від контакту розглядаються, як 
рівноправні, два вирази, один з яких, а саме, 


І ї І -иІХ х)- х)- 
у соучм М ОМе (2), Ас(х)- с(х)-в, (6) 


ху) з? 


являє собою лінеаризацію зв'язку 





що б 
- а х) -. 45 зиц(х), (7) 
пли. (6-х) КС Є) 
а другий - зв'язок вигляду (5). Так чи інакше, але в наступних роботах О.К. Маловичка | 14, 15), В.М. 
Новоселицького (16, 17) і багатьох інших авторів |7, 18-34) за рівняння для контактної поверхні брали 
вираз (7) або його лінеаризацію (6). Кожен з цих виразів залежить від двох параметрів - с1 й - Й описує, 
як це випливає з інтегрального рівняння (4) Заморєва при а-» -о0 і р -» со, вплив необмеженого в ла- 
теральному напрямку і обмеженого за товщиною геологічного тіла на зразок синклінорію, антикліналі й 
синкліналі якого виповнені масами однорідної густини протилежних знаків (рис. 3). Тіло це відділено 
від рівня спостережень шаром товщини 1, дія якого виразом (7) не враховується. Разом з тим, у цій пос- 
тановці було отримано цікаві результати, які бажано використати в теорії структурної задачі, що грун- 
тується на рівнянні (3). 
До числа таких розробок в першу чергу відно- 
сяться відкриті М.Г. Сербуленком |18, 191 і 
Ю.В. Антоновим |20-22)| способи визначення 
двох і більше контактних границь, що розділя- 
ють однорідні шари. Слід згадати також дос- 
лідження О.К. Маловичком (15) питань уточ- 
нення розв'язків лінеаризованого рівняння та 
єдиності його розв'язків, якими встановлено, 
що одному й тому ж полю можуть відповіда- 
ти як обмежені однорідні випуклі тіла, так і 
різноманітні необмежені однорідні тіла типу шару. Заслуговують на увагу перші результати чисельно- 
го розв'язку нелінійного інтегрального рівняння (7) у межах теорії регуляризації |23-301, вивчення В.І. 
Старостенком (31) можливостей стійкого чисельного визначення декількох контактних границь, що 
описуються нелінійними рівняннями виду (7), 1 спроби О.І. Кобрунова та Р.П. Денисюка |32-34) відшу- 
кати за допомогою одного чисельного алгоритму границі розподілу неоднорідних комплексів. 
Аналогічну модель задачі для контактної поверхні в комплексній площині 2-х іу за- 





Рис. 3. Модель Маловичка. 


пропоновано в 70-х р.р. минулого століття В.М. Страховим |35-39). 


4 
Саме для шаруватого середовища у смузі, 
відок ремленій від нуля і обмеженій прямими 


узі у2й", й «й «0, що складаються з 
двох однорідних шарів порід з густинами с, 1 62, 
які розділяє жорданова крива дОРД:пео п), 


П9СХОЮ зі стрибком Є - 2, - б, густини на 





контакті і з асимптотою й, й «й«й"', однією й 
Рис. 4. Модель Страхова. тією ж ліворуч та праворуч, яку крива ОЇ) може 

перетинати не більше, ніж скінченне число разів 
(рис. 4), він вивів нове нелінійне інтегральне рівняння 


у 


зобі ; р 
з 5 ПДР и" (2) іє (8) 
и ра и ( , ієр 
де Р)" - частина площини "вище" кривої С), )' - "нижче? кривої ДО), з3-2чін, 5-2-їй; ца) - 


в калібровці (3) задано за результатами спостережень на "земній поверхні" у - 0. Для цього рівняння 


В.М. Страхов довів теорему єдиності його розв'язку і редукував задачу визначення контакту за полем 
до задачі побудови допоміжної функції, що конформно (і однолистно) відображає нижню півплощину 


певної комплексної площини на область /)? під шуканою кривою. 

Дей підхід до розв'язування задачі передбачалося побудувати за аналогією з розвинутим В.К. Іва- 
новим методом відшукання однорідно розподілених джерел на скінченних однозв'язних областях |40, 
41). Незабаром після вибраної редукції було знайдено цілі класи розв'язків рівняння (8). Їх подавали або 
як сімейства розв'язків у скінченному вигляді, коли границю шуканої області визначають скінченним чи- 
слом параметрів зі скінченної системи (нелінійних) рівнянь, або у вигляді двопараметричних сімейств 
розв'язків "в малому", коли границю визначають у вигляді (невеликих) відхилень контура шуканої від 
контура заданої області. Ці класи отримали, відповідно, О. В. Цирульський |42, 43| 1 В.Г. Чередніченко 
144). До того ж В.Г. Чередніченко зняв обмеження |38| на розв'язки рівняння (8), допустивши, що крива 
ОР необов'язково має бути графіком і може перетинати асимптоту Й нескінченну кількість разів. 

Отож, для опису фундаментальної залежності гравітаційного поля від шаруватого середовища 
було запропоновано на початок 1980-х років дві принципово різні математичні моделі: модель Нумеро- 
ва (рис.1) та модель Маловичка-Страхова (рис.4). У моделі Нумерова як прикладну модель середовища 
розглянуто необмежений шар між контурами 0), та ОР,, його притягання враховано у вигляді суми 
притягання шаром постійної товщини / і відхиленнями від цього шару між прямою у-Л і кривою 
зроеній (х). У моделі Маловичка-Страхова предметну модель шаруватого середовища ототожнено з флук- 
туаціями у з л(х)- й кривої у - л(х) відносно своєї асимптоти й, що обмежують "додатні" та "від'єм- 
ні? маси густиною с; поле цих мас аж ніяк не могло збігатися з полем моделі Нумерова, оскільки не бу- 
ло враховано вплив шару товщиною Й. Проте кожне з інтегральних рівнянь (1) і (8), що відповідало тій 
чи іншій моделі, виявилось однаково залежним від двох параметрів: надлишкової густини 6 1 сталої то- 
вщини (асимптоти) / шару. Відповідно і розв'язки рівнянь також стали залежними від цих параметрів. 

У літературі відзначений параметричний зв'язок подано навіть як визначальний для класифікації 
обернених задач гравіметрії. З виходом у світ раць (35-39, 42-44), в яких закладено основи теорії обер- 
неної задачі для предметної моделі Маловичка-Страхова (рис.4), устоялося положення, що обернені за- 
дачі для обмежених тіл вринуцново відрізняються від обернених задач для необмежених тіл. Їх відмін- 
ність полягає саме в тому, що розв'язки задач для обмежених тіл залежать тільки від одного параметра 
а, тоді як розв'язки для необмежених тіл залежать як від параметра а, так і від параметра Й ("середньої 
глибини" за Нумеровим або асимптоти за Страховим). 

І довгий час ніхто не піддавав сумніву переконливість цього твердження, хоч воно явно не узгод- 
жувалося з аналізом рівнянь (8) та (5). Уже побіжний погляд на них переконує, що рівняння (8), а отже, 
й його розв'язок має залежати як від параметра є, так 1 від асимптоти й шуканого контура 0), тоді як 
рівняння (5) для того ж контура (як і численних інших обмежених) від параметра Й не залежить. Рівнян- 
ня (5), якщо воно вірне, давало зразок опису задачі для контактної поверхні, для якого одне з положень 
теорії, розвинутої в роботах (35-39, 42-44|, потребувало уточнення. З невідомих причин цього не стало- 
ся. Або автор пропозиції (5) чи його послідовники не були впевнені в адекватності опису моделі задачі, 
або авторам праць (35-39, 42-44) не була відома праця (131. 

Так чи інакше, але рівняння (3) було проігноровано під час побудови теорії обернених задач для 
шаруватих середовищ і питання про вирішальну відмінність задач для обмежених та необмежених тіл 
залишилося відкритим. Окрім того, в створеній теорії не розглянуто узагальнення відкритого Г. Рейнбо- 
єм факту, що розв'язки лінеаризованого інтегрального рівняння для контактної поверхні слід шукати на 


5 
обмежених множинах гармонічних функцій. Математичне обгрунтування згаданих та інших супутніх 
проблем стало предметом подальших досліджень |45-551. 

Початок цих досліджень було закладено аналізом предметних моделей Нумерова і Маловика- 


Страхова. В результаті адекватною для двовимірною евклідового простору ко) - В) х кі ; 
В) - (- со, со) виявилась предметна модель Нумерова, яку подано у вигляді шару 

р- (ху) -о«х«о; 04ук«4х)), 
обмеженого контурами ДІД): у-0 і ОД,: ус (х), хе Ек), зануреного в простір КО) і виповненого 
масами постійної густини с - о, - с» . Відповідну предметній математичну модель контактної поверхні 


знаходять за допомогою регуляризації відповідного означеного інтегралу, що розбігається на границі 
Ор «др, 00р, збурювальної області Д. Регуляризацію здійснено за правилом |45, 52| 


пр у ДРУ ЗСУ - 
НО 1 олені; урей 
-щх), ус0, 





І -х)з -уЇ 
Паб НЕЇ ше 200-2у, Оу «Сб 
- ху-у 
- (2), 2 С(х) 
Ця операція дає теоретичне обтрунтування нелінійного інтегрального рівняння (5) і забезпечує 
його адекватність геологічній моделі шаруватого середовища. Аналог цього рівняння у комплексній 
площині 2 - х -іу має такий вигляд |45, 53| 


55 1 НН ас и(2), (9) 


2 27. р, 5-7 





де ца) - напруження поля за калібровкою (3); 5-2 кій (2); зп ХНіС (х) - точки контура ОЙ),; 
5 є-4612), 5, ЗХ - іс (х) - спряжені з ними точки. 

Перетворимо це рівняння до зручного для порівняння з поданням (8) вигляду. Оскільки, як зго- 
дом буде доведено, 





і (10) 
гр, 





і аз чіл, іє р, 
іл, іє Ї), 
де Д" - область над кривою 0),,а Її) - під ОЇ, то, очевидно, 


М, Ї 5-5- з -5а) Зо) зц(2). 


длі і, 5-1 


Якщо ж серед множини обмежених знайдеться крива ОП, , що має асимптоту Й ліворуч й право- 


руч, то, взявши у співвідношенні (9) 5 - 5, З -ни іс (х) -2ій, отримаємо 
хо» 

и о ог зи(2) або ж «о Ї 

длід, | 8573 длі з, 
Порівнюючи останній вираз з рівнянням Страхова (8), дійдемо висновку, що рівняння (8) є частинним 
випадком від наведеного (9), в якому до того ж не враховується (як і у попередньому рівнянні (7)) при- 
тягання шаром постійної товщини 1. Окрім того, подання (8) від рівняння (9) суттєво відрізняється не 
тільки тим, що воно описує значно вужчу множину можливих контурів збурювальних шарів, а й, голов- 
не, тим, що воно залежить від параметра 1, якому в теорії |35-39, 42-44) надано визначальну роль. 


5 -5 - ЗіЙ 


5-72 


аз ца)-ній й 


У зв'язку з вивченням шаруватих структур введено |45, 47| клас Нумерова Ми (с П) Іо 
збе (р)х С (ка) б) як декартовий добуток множини 5), (р) обмежених і (майже всюди) локально ін- 
тегрованих функцій, що описують густину с « с(х, у), ба у)є Р мас в області Р с В? та сукупності 
сне (ві9) Кк -0,1,..,; 0«сх1 функцій у- 202), х є В? неперевних (за Гельдером) з їхніми похід- 
ними до К-го порядку включно з показниками о, що описують контакти різнорідних шарів зі смуги 
Па Со, у): -о«х«о; 0 СО) « т де П, - обмежена стала. У працях |45-55| розв'язки рівняння (9) 


є г - й 
розглянуто на класі Ми! «і, п) зі сталою густиною (відзначеною символом 1). Показано |45|, що на 


класі обмежених східчастих контурів обернене відображення поля в параметри середовища є розрив- 
ним. Більше того, оскільки |46| нелінійний оператор 








157, (2-х) нсх 
4сх)- Є00-з- | не; 
ях (8-х) чє"(з) 
прямої задачі для контактної поверхні є обмеженим, неперервним і компактним тоді 1 тільки тоді, коли 


(9 


він діє на множині Ми"(1, П) функцій, неперервних (з показником а) на К"" разом зі своїми першими 


похідними, то задача розв'язування рівняння (5) поставлена некоректно і на множині Мибе(1, п) непе- 


рервних функцій (цей факт, між іншим, випливає також з того, що кожна східчаста функція може бути 
подана у вигляді границі певної послідовності неперервних функцій). 

Доведено теорему єдиності розв'язку рівняння (5) на Ми (1, п). Доведення здійснено узагальне- 
ним на випадок необмежених зірчастих тіл методом, створеним зусиллями П.С. Новікова |56|, О. О. За- 
морєва (11), Л.М, Сретенського |357|, І.Т. Тодорова, ДП. Зідарова (581, М.М. Лаврентьєва (59) 1 О. І. 
Прилєпка |60| для обмежених зірчастих тіл. Суть методу полягає в побудові й оцінці певного функціо- 
налу, значення якого слугують ідентифікатором дійсності висловленої пропозиції. Функціонал будують 
з використанням фундаментальної інтегро-диференціальної тотожності Остроградського-Г ріна, він ві- 
дображає рівність зовнішніх потенціалів (або їх похідних) і властивості густини тіла, що не породжує 
зовнішнього потенціалу. Функціонал визначають на певній гармонічній функції, яку, у свою чергу, бу- 
дують за визначеними даними граничної задачі Діріхле в об'єднанні областей, що зумовлюють згадані 
потенціали. 

Вибором методу доведення теорем єдиності розв'язку визначають 1 сам структурний клас 
Ми (1, П) , до якого цей розв'язок може належати. 

Проблему же існування розв'язку задачі, як її трактують у працях |47, 52-55|, а також питання 
стійкості розв'язків розкривають майже "дарма" завдяки спеціальному поданню лівої частини рівнян- 
ня (5). Задавши її у вигляді двох складових, одна з яких (невідома фунція), як виявилося, домінує за нор- 
мою над іншою (інтегральним членом, залежним від невідомої функції), автори отримали можливість 
побудувати таку послідовність компактних операторів, яка збігається до необмеженого оберненого до 
ИТГа 59 оператора, 1 за Його допомогою, у свою чергу, визначити точний розв'язок рівняння з множини 


допустимих. 
Таким чином було встановлено, що задачу відновлення контуру необмеженого шару за даними 


5 : : 1 
спостережень потенціалу поставлено коректно за Тихоновим на нумерівському класі Ми! ЗП) , тоб- 
то, що на класі ми, П) розв'язок задачі існує, єдиний 1 стійкий, До того ж, цей розв'язок залежить 


(як 1 в задачах для обмежених тіл) тільки від одного параметра с і зовсім не залежить від параметра Й 
(середньої глибини до контакту за Нумеровим або його асимптоти за Страховим). 

На основі створеної теорії розроблено два чисельні регулярні способи розв'язування оберненої 
задачі за наявності перешкод. Побудова одного з них заснована на тому, що значення поля, за якими ви- 
значають фрагмент контуру, задають на практично необмеженому інтервалі. В цій ситуації апроксима- 
цію оператора прямої задачі АС ;х) можна виконувати з високою точністю. Побудову другого здійсню- 


ють за умови, що значення поля задано на короткому інтервалі і апроксимацію оператора прямої задачі 
можна здійснити з невисокою точністю. У першому випадку, за умови високоточної апроксимації опе- 


ратора ИТуя о. для чисельного пошуку каркасу наближеного розв'язку рівняння (5) конструюють регу- 


ляризуючий алгоритм за способом М.М. Лаврентьєва |591; у другому - для знаходження каркасу розв'я- 
зку використовують алгоритм, побудований в межах теорії регуляризації за Тихоновим (29). 

У цьому алгоритмі вперше обтрунтовано вибір саме такого стабілізатора, який відповідає хара- 
ктерним властивостям розв'язків й забезпечує пошук розв'язків на множинах єдиності, і рекомендовано 
певну нестандартну процедуру для точнішої порівняно з традиційними апроксимації рівнянь Ейлера, що 
не потребує перемноження матриць високої розмірності. Ретельно вивчені питання скінченновимірної 
апроксимації оператора УТГа ро) і послідовних наближень оберненого оператора Ач ;х). 


Особливу увага приділено залежності чисельних розв'язків від багатьох параметрів, які супро- 
воджують дискретизацію задачі, і вказано способи вибору таких значень цих параметрів, які сприяють 
одержанню оптимальних чисельних розв'язків. Доведено збіжність каркасів наближених розв'язків до 
точного розв'язку рівняння (5). Детально з теорією оберненої задачі для контактної поверхні читач мо- 
же при бажанні ознайомитись в дисертації Н.Н. Чорної (521. 

У працях (61-70) структурну задачу гравіметрії трактовано як задачу визначення за даними спос- 
тережень контуру необмеженого зірчастого тіла, близького до шару постійної товщини. Основним рів- 
нянням задачі в цій постановці є співвідношення (1), одержане Б.В. Нумеровим понад 70 років тому. 


Для відновлення границі О),: у б (о) хе В? збурювального тіла відомої сталої густини с за значе- 
ннями похідних логарифмічного потенціалу, заданими на множині єдиності, або, що те ж саме, для роз- 


ї/ 
в'язування нелінійного інтегрального рівняння (1) запропоновано |6Ї, 64) такий спосіб. Розв'язок рівня- 
ння (1) відшукують у вигляді границі послідовності розв'язків певної послідовності лінійних інтеграль- 
них рівнянь першого роду, що апроксимує вихідне нелінійне рівняння. Встановлено |62, 65|, що опера- 
тори послідовності лінійних рівнянь є обмеженими, неперервними і компактними в просторах функцій, 
які найчастіше використовуються при розв'язуванні обернених задач потенціалу. 

Виконано спектральний аналіз операторів, на основі якого вказано ефективні критерії належнос- 
ті вихідних даних 1 розв'язків відповідним Їм області значень оператора 1 області визначення нормаль- 
них розв'язків. Доведено |63, 66| локальні теореми існування та стійкості розв'язків задачі визначення 
необмежених зірчастих областей, близьких до заданих. Показано, що задачу поставлено коректно 
тільки на компактних множинах банахових просторів вихідних даних та розв'язків. Внаслідок умовної 
її коректності для отримання ефективного розв'язку слід застосовувати спеціальні регулярні методи. В 
основі підходу закладено теорію згладжувального функціоналу над просторами даних 1 розв'язків. 

Детально вивчено його складові: функціонали нев'язки або функціонали типу нев'язки та стабі- 
лізатори |67|. Показано, що функціонали нев'язки чи типу нев'язки лінійного інтегрального рівняння 
першого роду з невід'ємним компактним оператором є строго випуклими, двічі неперервно диференці- 
йованими функціями. Множини точок мінімумів функціоналів над Гільбертовим простором даних є 
множинами другої категорії по відношенню до множини всіх точних (нормальних) розв'язків інтеграль- 
ного рівняння (множини єдиності). Тому не кожна послідовність, яка мінімізує функціонали, збігається 
до множини єдиності, через що задача пошуку мінімуму функціоналів над гільбертовим простором да- 
них є некоректною. 

Недостатніми є й додаткові вимоги, які широко використовують, стосовно мінімальності норми 
розв'язку при мінімізації функціоналів нев'язки. Для арентування пошуку нормального розв'язку рів- 
няння за допомогою мінімізації функціоналів нев'язки й типу нев'язки запропоновано спеціальні стабі- 
лізатори, визначені на множині єдиності. Їх конструюють за допомогою певних диференціальних опера- 
торів, власні функції яких тісно пов'язані з власними функціями операторів прямої відповідності в зада- 
чах визначення зірчастих областей. Спеціальні стабілізатори - - це невід'ємні, сильно випуклі, двічі не- 
перервно диференційовані функціонали, що не мають локальних мінімумів, окрім глобального. 

Вперше встановлено існування такого регуляризуючого оператора для лінійного інтегрального 
рівняння першого роду, який серед відомих його модифікацій має найкращу обумовленість. Кожній 
функції з гільбертового простору даних оператор ставить у відповідність функції гільбертового просто- 
ру розв'язків, що мають узагальнену похідну, яка інтегрується з квадратом. У свою чергу, кожна з функ- 
цій, яку обчислюють за допомогою цього оператора, є екстремаллю згладжуювального функціоналу, 
який є зваженою сумою функціоналу типу нев'язки з одним із запропонованих стабілізаторів, що забез- 
печують однозначну розв'язуваність вихідних задач. Глобальні регуляризуючі оператори визначення 
зірчастих областей, близьких до шару постійної товщини, будують у вигляді послідовності локальних 
регуляризуючих операторів для відповідних послідовностей лінійних інтегральних рівнянь першого ро- 
ду, в правих частинах яких використовують згладжені вихідні дані. У свою чергу, локальні регуляризу- 
ючі оператори для лінійних інтегральних інтегральних рівнянь конструюють у вигляді параметризова- 
них послідовностей, що мінімізують згладжувальні функціонали |б6381. 

Ця стаття розпочинає серію, в якій передбачається розгорнути теорію структурних задач, яку бу- 
ло започатковано в нотатках |45, 48|, зокрема розглянути питання коректного розв'язування рівняння 
(9) в комплексній площині. 

Проблеми коректної розв'язуваності задачі в комплексній області. Розгортання теорії роз- 
в'язування рівняння (9) розпочнемо з обгрунтування співвідношення (10). 


Лема. Якщо контур ДО: у-б (х), х є В) належить до класу неперервно диференційованих (за 
Гельдером) функцій Ми(1, П), то 





1 | аз 91, є р, 
-І -1, сіє р 
Доведення проведемо для випадку, коли 2 є ДЮ". Оскільки функція виду  / (з) 1 (з -2) зі 


аналітичною в незамкненій області ПД" Ма) з виколотою точкою 2-х кіу, то на основі відомої теоре- 
ми Коші |71| інтеграл, що нас цікавить, можна подати в такому вигляді (рис. 5): 


аз аз 1. аз 
а вай Ці Пе ре -І | РОН М 


ха ЗОМ Од др» хьузчоМ С І 2 з са тС8 


де др) ід); -- частини контуру ОО між точками (х,, Сх, )), (0,2(0)) ї (2. Сх, )); С, - розтин між 
точками 2 з (0,0) і (0, 2(0)); С, і С. - коло радіуса є - і і 























5-2 
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В- Я 2 з центром у точці 2 - (0,0); Сі С, - дуги 


секторів, кути розхилу яких відповідно дорівнюють: 














Ф, 7 агсіє сх) - з агсіє Щ - »0, 
| ГРА НН 

Ф, - агсіє сх) шо Р тонна ннннно -»0. 
Ж 1з- (о, / 2 (х, )) ре 








(11) 


б. схил ую ж . 
Криволінійні інтеграли за дугами ОЇ), і- 12, 


І 





очевидно, існують тільки тоді, коли контур 0 (х) нале- 


у 50) 
: п : 1, 
Рис. 5. Зображення інтегралу для нелінійного рів- 2жить до класу Ми! «(1 п). У цьому випадку величину 
няння контакту в комплексній площині. 








1 аз о 1р о 1чіє" Чо зніє" 
Ї сл . Ї ІЗ ) ає | і") 42 
грічері 9 Лі мода 43) а. печіС(Є)- 2 

цілком визначають фіксованими числами х,, і х,, а її границю при х, -2-00, х, -» 00 з урахуванням 


відношень (111) знаходять так: 


2. їф л .- і 2л -їФ ляФ, .їФ 
1 й і 4ф 2. 4 | Кіе 4 | Кіе 4 М 
З є"? усне зм К-е" К-е" 





Випадок, коли , розглядають аналогічно. Лему доведено. 
На основі доведеної леми нелінійне інтегральне рівняння (9) подамо у вигляді 


(5 зно Гм) (12) 


ор 





зручнішому для подальшого дослідження. 
Теорема І. Якщо контур Ор: у Сх), хе В належить до класу Ми, П), тобто 








"ЩЕ М «т, т, М «со, то оператор прямої задачі 


Мега чно | 


опаераоа 


задовольняє вимогам | (о з 


5-5 


аз 





. 1 
є обмеженим, неперервним і компактним оператором на С (кі 7) 
Доведення обмеженості оператора випливає з очевидної оцінки 


УУЗВ В РУКОЮ ССЗ РА РЕТРО 
до У По5-2 


ор 


«3М. 














Для доведення його неперервності розглянемо різницю 


1 с5-5 5 1 г8-5 
Ас трг)- Аг) є «діб Р а Га, 


45 783 2лі др 


де 5 збчідчіпаячіт, др": уз 2Сх)ч л(х), хе В); пох) -- достатньо мала варіація функції С) 





у тому сенсі, що функція п(х) для довільного числа є » 0) задовольняє вимогам 
|У, - тлахриоо) «є; (-о2)- т/оо)- 0. 


Оскільки варіація п(х) мала, то 


і кока 
1 ад 1 Де С "Е8 З 
5-2 5-21чіц(в-1) 5-7 5-2 158-7 








і тому АС з та)- А(2 2) М с Ї ей а я я Я С ітріз Р. 
Ул Дн од 2лі 55 - 1) 
Звідси неважко вивести рівність (С з щі)- Й (а :2) 2 Со «Т7є/2, яка засвідчує неперервність 


2 


оператора ИТГа ,2) вс (ві ) через те, що число є » () довільне. 
Для доведення компактності (повної неперервності) оператора ИТ/а 52) покажемо, що він 
переводить довільну обмежену множину, зокрема клас МиеУІ, П) , функцій в компактну (за метрикою 


С (в?) множину функцій ца) АС :2)є У(А). Раніше вже була нагода переконатися, що оператор 


ИТГа :2) обмежений на довільній обмеженій множині функцій, тому для сє Ми, П) матимемо 


ме, «зіва 
Цією нерівністю встановлено рівномірну обмеженість множини У(4). Доведемо її одностайну 


неперервність. Нехай Їх, я 25 « є для довільного числа є » 0. Тоді різниці 


Е 1 ну б 
ца) - ма, ) є 2ЙС Со) вх - ПРЕ и 2 
з урахуванням розгортання 
п пі 
1 - 1 1 1 раль, асо, І 2-8, 
дод. чо, Ї (а 2,)/(5 2,) 8чзіві зів, зоб, о бот, 
можна надати вигляду 


ца) аз) о (а) ооо, зо) сій ко) (2 -2,)" аа М 


2лі раї 53 - 2,) 23 - 2 У» її 2, у" 


Проте функція 0 (х) неперервно диференційованна на Ек і, до того ж, |Х, - х, « о Й 





І 
а 

















, Через це норма 








різниці набирає вигляду: 


Іа) -м(2. | Їм є і ; , Т, 





5-5 
з ав| « є(дт- 1/2) 
то 
Одержана нерівність означає, що множина И(А) функцій ца) є одностайно неперервною. За ві- 
домою теоремою Арцела |72|, множина У(А) є компактною, а оператор 42 ,2) - цілком неперервним 


(компактним) у метриці простору неперервних функцій. Теорему повністю доведено. 

Визначальна роль доведеної теореми для розв'язування рівняння (12) полягає в тому, що вона 
орієнтує напрям пошуку розв'язків за їх характеризацією. На конкретних елементах того чи іншого кла- 
су функцій можна показати, що неперервна залежність розв'язку інтегрального рівняння (12) від його 
правої частини відсутня не тільки в просторах обмежених (східчастих) функцій, а й у просторах непере- 


рвних С (в) або ж інтегрованих (абсолютно чи з квадратом) І9(к9) Кк - 1,2 функцій. 

Для того, щоб у цьому переконатися, досить розглянути значення оператора АС ,2) на певній по- 
сліджовності неперервних функцій, що збігаються до тієї східчастої, яку розглянуто у праці |33|. Опера- 
тор АС 2) , Як встановлено теоремою 1, неперервний принаймні на класі Ми? (1, П) , і щонайменше ті- 
льки на цьому класі потрібно відшукувати розв'язок рівняння (12) для вхідних даних з просторів С (в) 
та Ук?) к12. 

Теорема 2 (єдиності). Розв'язок нелінійного інтегрального рівняння (12), якщо він існує, єдиний 
на множині Ми) (п) неперервно диференційованих функцій. 
Дійсно, однорідне нелінійне рівняння 











ПЕРЛИ ЗОРЕ рН (13) 
поь 
на множині ми (1, П) має тільки тривіальний розв'язок. У цьому нас переконує суперечлива нерів- 
ність 
1 г5-5 1 1 го аз 1 
-уЇ з ах « 5 Ї « я, 
нм Бо Пот кре ГЕ ер, 


яку виведено на основі леми. 
тже, однорідне рівняння має тільки тривіальний розв'язок, а рівняння , якщо Й має 
о 13 ' 12 
розв'язок, то він однозначно визначається на класі Мич1, п). 


Надалі зауважимо, що з огляду на те, що довільну функцію 0 (х) з множини Ми (1, П) можна 
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подати у вигляді 2 (х)- С(0)-- Їс (г)4т, ця множина є рівномірно обмеженою і водночас неперервною, 
тобто є компактною знана функцій. 

Теорема 3 (їснування). Процес послідовних наближень 
де і доні й А. і 5-х) 


2лі аб» 


п 0 -(о 
рус занен 45" - ца), п є 0,2,уг 8-5 аг) 4) 
) 
збігається на компактній множині ми, П) до (нормального) розв'язку нелінійного інтегрального 


рівняння (12) зі швидкістю геометричної прогресії. 
Для доведення цієї теореми достатньо показати збіжність ітераційного процесу (14) і визначити 


його швидкість. Твердження стосовно того, що гранична функція 0 (х) послідовності Я. о) буде 


розв'язком нелінійного інтегрального рівняння (12), випливатиме з теореми єдиності розв'язку. 
Починаючи з'ясування збіжності послідовних наближень, розглянемо різницю 











ло «ок п чт 1 50?) го 5"? п 1 517) ри 5141) п- 
(зі ву) з ее | пут Злі | су дхтг) 
РАБ ОТ о РН м 
Ршясспрошення: икладок знозначимо зон) о) Це 0) ро), п -- 0,1,2,..., 1, скориставшись 


тим, що прирости по) за метрикою простору неперервних функцій С (кб) набагато менші самих наб- 


лижень 0 (п) контуру 5, матимемо 


ГОД І сне ЛІХ са р Вл 
зі; 50 а 1знітіО Д5)- 2 50Ї-г зай МО Р 


ахтЯ) що дз іані?), п - 0,1,2,... 
У цих позначеннях попередня рівність виглядає так: 


(з) (я) (а) и) 
нова и і 5 біо і ен дн | ні РОТИ 











пий о ор Злі дн) 5002 Злі 5-7 
(а-1) (а-1) 
Я 1 б дах? се Ї я ах'ттр 
Злі (55-38 Злі 1.05 














(п-1) (п-1) -(п-1) 
у. 1 2 (ау 8 8 ані) 
і - | Р дере 4, (з нії г) | ; 


І Ж 2 (п-1)  «п-і) Ж 
еЗ з "Ю РЕ ще се я у 


Ця нерівність слугує основою для наступного ланцюжка 
а от і. « 229 зані « ЄМ «У «мб ай є 


який, у свою чергу, засвідчує, що оператор АС 52) задачі для контактної поверхні є стискаючим. Наб- 


з якого встановлюємо нерівність 


(р 











, 


грі 


«ре 























ей 


лиження х) п 0.1,2,..., збігаються до розв'язк х) рівняння зі швидкістю геометричної 
с), п 0,2,» зб розв'язку Сх) р (12) р 


баб . чай 
прогресії зі знаменником 2 ". Теорему доведено. 
Зауваження. Оператори 
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послідовності Ви, С ка), яка, очевидно, апроксимує при 500) -з 510) ш ца) обернений (нелінійний) опе- 
ратор Р я 12) рівняння для контактної поверхні, є обмеженими, неперервними 1 компактними опера- 
торами на компакті И(А). Виявити ці властивості можна за допомогою засобів, які було застосовано під 
час досліджень оператора ИТГа 52) прямої задачі. До того ж, кожен з операторів Ви, С (9); 2) побіЗ 
проектує на нумерівську множину Ми а (1, п). Дійсно, для переконання в істинності останнього 
висловлювання досить у зображенні оператора перейти до параметричного його подання, а саме 





а ббдлні, 
Ви, Ся з о Є -и(2 
«1 5 ніс"(5)-2 і 
п) 
Наявність у зображенні (інтегрованої) похідної аа від наближення 0 є) контуру 0 (2) МР а Ек), 


є ознакою функції, що належить до множини миле, п). 

Отже, в перебігу доведення теореми існування розв'язку знайдено таку послідовність 
Ви, С а) компактних операторів Ви, (2) пе 0,1,2,..., яка обмежену множину (А) функцій ца) 
однозначно проектує на компакт Ми, П) і прямує за метрикою простору неперервних функцій 
ско) на певному елементі ца) множини У(А) до нелінійного оберненого оператора АЧС :2) 
рівняння для контактної поверхні. 


Теорема 4 (стійкості). Нехай функції С Ох), хє В), і-12, що описують границі розподілу 


і 
однорідних шаруватих середовищ, належать до класу Мио(1, п). Якщо за достатньо малого числа 
є» 0 поля иа), Ж- Ек, які породжуються відповідно контурами оС, (х), підкоряються вимозі 

Ім (2)- а (2, «є, 
тобто їх вважають близькими, то самі границі також незначно відрізняються одна від одної в сенсі 
виконання умови | і (х)- із 9 Р-ЗУ/оя 


Для доведення теореми розглянемо різницю 


(2) - 4. (х)є з 39 аз? - Ї с) аз? ш(2)чн (2). 














дл а 59-; 
Оскільки 

1 2(2)-65(2) (у, 1 3, (у 1 с3З| о 

т, кора а че еаЙлове з 
1 р 8-8) до Ї о УНК пнісо) аЄ «(9-62 х 
лід 87І Теніс )-2 ні ()-х || дя 9 

аз) Ра Р 459 450) 1 
тая оз еввонх 5100-65 





то Ра (х)- 24 ТЕЗИ з З (2)- и, (2 «2є ітеорему доведено. 


Доведенням теореми стійкості завершено дослідження умов коректної розв'язності оберненої за- 
дачі логарифмічного потенціалу для контактної поверхні в комплексній площині. Задачу поставлено ко- 


ректно на компактній множині МиеУ(І, П) неперервних разом з першими похідними функцій саме то- 


му, що на класі Ми (п) розв'язок задачі не тільки існує, а водночас є є диним і стійким. 

Особливістю проведеного теоретичного аналізу задачі є редукція задачі до розв'язування неліній- 
ного інтегрального рівняння (12) типу рівняння Урисона. Це рівняння отримано відповідною корекцією 
(на основі визначальної леми) рівняння для контактної поверхні в комплексній площині, виведеного у 
статті (45). Ліву частину отриманого рівняння, що залежить від невідомої функції, подано у вигляді 
двох складових, одна з яких - невідома функція - домінує (за нормою) над іншою, що є інтегральним 
оператором над невідомою фунцією. Вигода такого подання (на зразок лінійного інтегрального рівнян- 
ня другого роду) чітко вбачається в напрямі дослідження рівняння. По-перше, подання виявилося досту- 
пним для вивчяення умов єдиності, існування та стійкості розв'язків винятково елементарними засоба- 
ми. По-друге, виглядом рівняння визначився і спосіб його точного розв'язання. 
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А з огляду на те, що точний розв'язок відшукують за допомогою певного процесу послідовних 

наближень, сконструйованого, у свою чергу на основі вихідного подання рівняння, це подання стає ви- 

рішальним і в практичній реалізації розв'язку задачі в умовах, коли праву частину рівняння задано з пе- 
решкодами на дискретній множині точок. 
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